Effects of the f/-boson on the inner edge of neutron star crusts 
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We explore effects of the hght vector [/-boson, which is weakly coupled to nucleons, on the 
transition density pt and pressure Pt at the inner edge separating the liquid core from the solid 
crust of neutron stars. Three methods, i.e., the thermodynamical approach, the curvature matrix 
approach and the Vlasov equation approach are used to determine the transition density pt with the 
Skyrme effective nucleon-nucleon interactions. We find that the pt and Pt depend on not only the 
ratio of coupling strength to mass squared of the [/-boson j p? but also its mass p due to the finite 
range interaction from the [/-boson exchange. In particular, our results indicate that the pt and Pt 
are sensitive to both g'^ / p^ and p if the [/-boson mass p is larger than about 2 MeV. Furthermore, 
we show that both g^ / p^ and p can have significant infiuence on the mass-radius relation and the 
crustal fraction of total moment of inertia of neutron stars. In addition, we study the exchange 
term contribution of the [/-boson based on the density matrix expansion method, and demonstrate 
that the exchange term effects on the nuclear matter equation of state as well as the pt and Pt are 
generally negligible. 

PACS numbers: 26.60.-c, 14.70.Pw, ll.lO.Kk, 97.60.Jd, 



I. INTRODUCTION 

The possible existence of a neutral weakly coupled 
light spin-1 gauge [/-boson which is originated from 
supersymmetric extensions of the standard model with 
an extra U{1) symmetry, has recently attracted much 
attention due to its multifaceted influences in particle 
physics, nuclear physics, astrophysics and cosmology. For 
instance, the [7-boson can provide annihilation of light 
dark matter which can be responsible for the excess flux 
of 511 keV photons coming from the central region of our 
Galaxy observed by the SPI/INTEGRAL satellite 
It is also proposed that the [/-boson can be mediator 
of the putative "fifth force" providing a possible mech- 
anism for non- Newtonian gravity, i.e., the violation of 
the inverse-square-law (ISL) of Newtonian gravitational 
force at short distance |6l-l23l| . Thus far various upper lim- 
its on the deviation from the ISL have been put forward 
down to fcmtometcr range [l^ [l^ [l^. Furthermore, 
the [/-boson can involve a rich phenomenology in parti- 
cle physics and nuclear physics and may have observable 
effects in particle decara [23-[27| and nucleon scattering 
processes |l6l . [isl Hol . [28l [29j , which can also put limits 
on the [7-boson properties. Studying properties of the U- 
boson is thus important for understanding the relevant 
new physics beyond the standard model. 

Very recently, the effects of the U -boson on the nuclear 
matter equation of state (EOS) and neutron star struc- 
ture have been investigated |30l - [33| and it is shown that 
the vector [/-boson can significantly stiffen the nuclear 
matter EOS and thus enhance drastically the maximum 
mass of neutron stars. In particular, by considering the 
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[/-boson, the stability and observed global properties of 
neutron stars can be reasonably explained by using the 
neutron-rich matter EOS with a supersoft nuclear sym- 
metry energy at supersaturation densities consistent with 
the available terrestrial laboratory data on the 7r~/7r+ 
radio in relativistic heavy-ion collisions from FOPI/GSI 
[s^ . [35| . while the supersoft nuclear symmetry energy 
at supersaturation densities generally can not support a 
canonical mass {1AM(^) neutron star if the [/-boson is 
not introduced [3l[. The [/-boson has also been intro- 
duced to describe the recently discovered new holder of 
neutron star maximum mass of 1.97±0.04 Mq from PSR 
J1614-2230 [iil using soft EOS's consistent with exist- 
ing terrestrial nuclear laboratory experiments for hybrid 
neutron stars containing a quark core described by MIT 
bag model using reasonable parameters [H, [131 , and it is 
found that the constraints on the U -boson properties are 
consistent with existing constraints from neutron-proton 
and neutron-lead scatterings [Tsl. [l9| as well as the spec- 
troscopy of antiproton atoms ,16*. 

In the studies about the U -boson influences on neutron 
star structure [sol - lssj , the exchange term contribution of 
the [/-boson to the nuclear matter EOS has been ne- 
glected and only the direct term contribution has been 
considered, leading to that the nuclear matter EOS de- 
pends only on the ratio the coupling strength to mass 
squared of the [/-boson, namely, / p^ ■ Physically, the 
exchange term contribution of the U -boson to the nuclear 
matter EOS will depend on both the coupling constant g 
and the U -boson mass p due to the finite-range interac- 
tion mediated by the U -boson. It is thus interesting to see 
how the exchange term contribution of the U -boson will 
influence the nuclear matter EOS. Furthermore, neutron 
stars are expected to have a solid inner crust surrounding 
a liquid core. Knowledge on properties of the crust plays 
an important role in understanding many astrophysical 
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observations [38l - l48j . The inner crust spans the region 
from the neutron drip-out point to the inner edge sepa- 
rating the sohd crust from the homogeneous hquid core. 
While the neutron drip-out density pout is relatively well 
determined to be about 4 x 10^^ g/cm'^ [i^, the transi- 
tion density pt at the inner edge is still largely uncertain 
mainly because of our very limited knowledge on the EOS 
of neutron-rich nucleonic matter, especially the density 
dependence of the symmetry energy [42l |43| . The transi- 
tion density pt and the corresponding pressure Pt at the 
inner edge might be measurable indirectly from observa- 
tions of pulsar glitches [H, |4^ . Since the J7-boson can 
have significant influence on the nuclear matter EOS, it 
is therefore very interesting to see how the J7-boson will 
affect the inner edge of neutron star crusts. 

In the present work, we investigate effects of the light 
vector J7-boson on the transition density pt and pressure 
Pt at the inner edge of neutron stars crust. The den- 
sity matrix expansion (DME) approach [HO, [HI is used 
to describe the exchange term contribution of the finite- 
range interaction due to the [/-boson exchange. Based 
on the Skyrme effective nucleon-nucleon interactions, we 
use three methods, i.e., the thermodynamical approach, 
the curvature matrix approach and the Vlasov equation 
approach to determine the transition density pt- As ex- 
pected, our results indicate that the pt and Pt depend on 
not only the ratio of coupling strength to mass squared 
of the [/-boson j p^ but also its mass p due to the finite 
range interaction from the [/-boson exchange. Further- 
more, we find that the pt and Pt are sensitive to both 

I p^ and p if the [/-boson mass p is larger than about 
2 MeV and both / p^ and p can have significant in- 
fluence on the mass-radius relation and the crustal frac- 
tion of total moment of inertia of neutron stars. We also 
demonstrate that the exchange term has minor influence 
on the nuclear matter EOS as well as the pt and Pt for 
the parameter values of g'^ / p^ and p considered in this 
work. 



II. THEORETICAL MODELS AND METHODS 

A. Nuclear matter symmetry energy 

The EOS of isospin asymmetric nuclear matter, defined 
by its binding energy per nucleon, can be expanded to 
2nd-order in isospin asymmetry 5 as 



(1) 



where p = p„ -I- pp is the baryon density with p„ and 
Pp denoting the neutron and proton densities, respec- 
tively; 5 = (p„ — Pp)/{pp + Pn) is the isospin asymmetry; 
Eq{p) = E{p,d = 0) is the binding energy per nucleon 
in symmetric nuclear matter, and the nuclear symmetry 
energy is expressed as 



In Eq. ([T]), the absence of odd-order terms in S is due to 
the exchange symmetry between protons and neutrons in 
nuclear matter when one neglects the Coulomb interac- 
tion and assumes the charge symmetry of nuclear forces. 
The higher-order terms in S are negligible, leading to 
the well-known empirical parabolic law for the EOS of 
asymmetric nuclear matter, which has been verified by 
all many-body theories to date, at least for densities up 
to moderate values (See, e.g., Ref. [12). As a good 
approximation, the density-dependent symmetry energy 
Esymip) can thus be extracted from -Esym(p) ~ E{p,^ = 
l)^E{p,6 = 0). 

Around the nuclear matter saturation density po, the 
nuclear symmetry energy i?sym(p) can be expanded as 



Esymip) — Esym{pn) + -ri- —] 

-J Po 



o{{p-^n (3) 

Po 



where L is the slope parameter of the nuclear symmetry 
energy at po, i.e.. 



L = 3p( 
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The slope parameter L characterizes the density depen- 
dence of the nuclear symmetry energy around normal 
nuclear matter density, and thus carry important infor- 
mation on the properties of nuclear symmetry energy at 
both high and low densities. 

The EOS of isospin asymmetric nuclear matter is a 
basic ingredient to determine the properties of neutron 
stars. For symmetric nuclear matter with equal fractions 
of neutrons and protons, its EOS Eo{p) is relatively well- 
determined. In particular, the incompressibility Kq of 
symmetric nuclear matter at its saturation density po 
has been determined to be 240 ± 20 MeV from analyses 
of the nuclear giant monopole resonances (GMR) |53l - 
\6^ . and its EOS at densities of 2po < p < 5po has also 
been constrained by measurements of collective flows (63j 
and subthreshold kaon production [65} in relativistic 
nucleus-nucleus collisions. On the other hand, the EOS 
of asymmetric nuclear matter, especially the density de- 
pendence of the nuclear symmetry energy, is largely un- 
known. Although the nuclear symmetry energy at po is 
known to be around 30 MeV from the empirical liquid- 
drop mass formula [66l. |67|. its values at other densities, 
especially at supra-saturation densities, are poorly known 
[53, [11]. During the last decade, significant progress has 
been made both experimentally and theoretically on con- 
straining the behavior of the symmetry energy at sub- 
saturation density and the value of L constrained from 
different experimental data or methods has become con- 
sistently convergent to about 60 ± 30 MeV [69l47lt (See, 
e.g., Refs. 73] for recent summary). Furthermore, 
the IBUU04 transport model analysis of the FOPI data 
on the 7T~ /tt"*" ratio in central heavy-ion collisions at 
SIS/GSI [U energies suggests a very soft symmetry en- 
ergy at the supersaturation densities [HI (See also Refs. 
[ZJ, [zH ) • These studies have significantly improved our 
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understanding for the EOS of asymmetric nuclear mat- 
ter, which have important implications on the neutron 
star physics. 



B. Skyrme-Hartree-Fock approach 

For the nuclear effective interaction, we use in the 
present work the so-called standard Skyrme force (see, 
e.g., Ref. [Z^), which has been shown to be very suc- 
cessful in describing the structure of finite nuclei. In the 
standard Skyrme Hartree-Fock (SHF) approach, the nu- 
clear effective interaction is taken to have a zero-range, 
density- and momentum-dependent form [76l |. i.e.. 



T~Lcoui is the Coulomb term given by 



Fi2(R,r) = to{l + xoPa)5{v) 

+ ii3(l + a;3^<xK(R)'5(r) 
6 
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+ ^h{l + xiP„){KH{v) + 5{y)K^) 

+ iWQ{ai+a2) ■ [K X J(r)K], (5) 

with r = r 1 — r2 and R = (ri -|- r2)/2. In the above 
expression, the relative momentum operators K = (Vi — 
V2)/2i and K' = — (Vi— V2)/2i act on the wave function 
on the right and left, respectively. The quantities P„ and 
CTi denote, respectively, the spin exchange operator and 
Pauli spin matrices. The a, to — t^, xq — are the 9 
Skyrme interaction parameters and Wq is the spin-orbit 
coupling constant. 

Within the standard form (see Eq. ([S])), the total en- 
ergy of the nuclear system can be written as: 



E = / 7^(r)dV, 



(6) 



with y, the Skyrme energy density. In the standard SHE 
model, the total energy density of a spin-saturated nu- 
clear system considered in this work is written as [76j 



n = IC + Ho+n3+H 



n 



fin 



'Hcoul 



(7) 



where IC = is the kinetic-energy term and Hq, "Hs, 
"He//, y-fin are given by 

Ho = to[{2 + xo)p^ - {2x0 + 1){pI + pl)]/^ (8) 
= i3p"[(2 + X3)p2-(2a;3 + l)(p^ + p^)]/24 (9) 
Ueff = [t2{2x2 + 1) - <i(2a;i + l)](r„p„ + TpPp)/8 

+ [ti{2 + Xi)+t2i2 + X2)]Tp/8 (10) 

■Hf^n = [3ti(2 + a;i)-<2(2 + a;2)](Vp)V32 
~[3ti{2xi + l) + t2{2x2 + 1)] 
X [(Vp„)2 + (Vpp)2] /32 (11) 

in terms of the 9 Skyrme interaction parameters a, to — t^, 
xo ~ x^. In the above equations, pi and are, respec- 
tively, the local nucleon number and kinetic energy den- 
sities, whereas p and t are corresponding total densities. 



(12) 

The nucleon single-particle energy can be obtained 
from minimizing the total energy of the nuclear system 
with respect to its wave function as 

., = |J^ + C/,= 2^ + C/;,. = n,p, (13) 

where Uq is the single particle potential while m* and 
U* represent, respectively, the nucleon effective mass and 
effective single particle potential, which can be expressed, 
respectively, as 



^2 1 



2mJ 2mq 



+ -p[h{2 + Xi)+t2{2 + X2)] 



+ -p,[t2(2x2 + l)-ti(2a;i + l)], (14) 



and 



= -to[{2 + Xo)p^{2xo + l)pg] 

+ ^<^hp''-'[i2 + X3)p' - (2x3 + 1)(P^ + pI)] 

+ ^t3P"[(2 + X3)p-(2x3 + l)p,] 
+ ^[h{2 + Xi) + t2{2 + X2)]t 
+ i[t2(2x2 + l)-ii(2a;i + l)]T, 

o 

+ ^[t2{2 + X2)-3ti(2 + xi)]W^p 
lb 

+ ^[3h{2xi + l)+t2{2x2 + l)]S/^Pq (15) 
16 

For protons, the additional Coulomb potential is given 
by 



1/3 



(16) 



For the kinetic energy density, we use in this work 
the results from the extended Thomas-Fermi approxima- 
tion [73, i.e.. 



rq^ap'/' + b^^^ + cV'pq, 
Pq 



(17) 



where a = |(37r2)2/3, b = 1/36 and c = 1/3. From the 
nucleon single-particle energy, the nucleon chemical po- 
tential in infinite nuclear matter, i.e., jj,q, can be obtained 
as the value of the single-particle energy without gradient 
terms at the Fermi surface — h{3'K'^ pqY^^ . 
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C. The weakly coupled light vector [/-boson 



The integrations in Eq. (|23|) are defined as 



Fujii [6| first proposed that the non-Newtonian gravity 
can be described by adding a Yukawa term to the con- 
ventional gravitational potential between two objects of 
mass nil and m2, i.e., 



Grnrrn2 

Vgra{r) = (1 



-r/A 



(18) 



with 



hq = / drr'p'sdK^) 



q Ji/Vq 



(24) 
(25) 



where a is a dimcnsionlcss strength parameter, A is the 
length scale and G is the gravitational constant. In the 
boson exchange picture, the light and weakly coupled 
vector [/-boson is a favorite candidate mediating the ex- 
tra interaction for the non- Newtonian gravity [l| , leading 
to the finite-range Yukawa potential between two nucle- 
ons which can be expressed as 



^UB(r) 



9 

9 e 



■fir 



(19) 



where g and p represent the C/-boson-nucleon coupling 
constant and the [/-boson mass, respectively. Compar- 
ing Eq. ([T^ with the Yukawa term in Eq. ([T5| . one can 
find the relations a ~ — g^/(47rGm^) and A = 1/fi (in 
natural units) where m is the nuclcon mass. By adding 
the Yukawa potential of Eq. ([Tn)) to the standard Skyrme 
effective nucleon-nucleon interaction in Eq. ([5|) , the extra 
binding energy of the nuclear system due to the U -boson 
can be expressed as the integral of energy density Htjb 
as follows 



E 



UB 



•HuB(r)dV, 



with 



Hub = ^uB + 



UB' 



(20) 



(21) 



where "Hub ^'^^ ^ub direct and exchange contri- 

bution to the energy density, respectively. For the finite- 
range Yukawa interaction of Eq. (|T9l) . the direct term 
contribution to the energy density can be easily obtained 

as fsO-lal 



n^B-^J Pin) 



g2 ^-tir 

47r r 



p{r2)dridr2 



\^^. (22) 



where V is the normalization volume, p = pn + Pp is the 
baryon number density, and r = jri — r2|. 

Although the direct term contribution of a finite-range 
interaction to the nuclear energy density can be treated 
exactly, it is numerically challenging to evaluate the ex- 
change contribution. The latter can be, however, approx- 
imated by that from a Skyrme-like zero-r ang e interaction 
using the density-matrix expansion [s^, l5l| and the re- 
sults can be obtained as 



«UB - 2^ T 



q—p,n 



2pqTgIlg + i 



dl 



lq+ Pq 



dpq 



i^Pq)' 



q—p,n 



li3n'r^'prhq+Pp2q 



(23) 



Pshikqr) 



ki 



-.hikqr), 



2(^^3(fcfr), 



(26) 
(27) 



where ji and are, respectively, the first- and third- 
order spherical Bessel functions and k^ ~ (Stt^p^)^/'^ is 
the Fermi momentum. 

One can see from Eq. that for the direct term, the 
[/-boson contributes to the nuclear energy density only 
through the combination g^/p^. On the other hand, it 
is indicated from Eq. that the exchange term con- 
tribution to the nuclear energy density depends on both 
the coupling constant g and the mass p in a, complicated 
way. Furthermore, the density gradient terms appear 
automatically in the exchange term contribution to the 
nuclear energy density due to the density matrix expan- 
sion of the finite-range Yukawa potential. It is interesting 
to see that the exchange term contribution to the nuclear 
matter EOS further depends on the isospin asymmetry 
although the direct term contribution to the nuclear mat- 
ter EOS is isospin independent due to the fact that the 
[/-boson is an isoscalar boson. Therefore, the [/-boson 
will contribute to the nuclear symmetry energy through 
the exchange term. However, as will be shown later, the 
exchange term contribution to the nuclear matter EOS 
and the symmetry energy is quite small and can be ne- 
glected safely within the parameter value region of the 
coupling constant g and the [/-boson mass p considered 
in the present work. 

The nucleon single-particle energy due to the U -boson 
can be obtained from variation of the corresponding en- 
ergy of Eq. ((20l) with respect to its wave function. The 
[/-boson contribution to the nucleon effective mass can 
be expressed as 



2™5,UB 



9 r 



(28) 



which should be added to the right hand side of Eq. 
to obtain the total nuclcon effective mass m^. It should 
be noted that the [/-boson contribution to the nucleon 
effective mass is from the exchange term with the density 
matrix expansion of the finite-range Yukawa potential, 
which leads to the isospin dependent contribution to the 
nucleon effective mass although the [/-boson is isoscalar. 
The [/-boson contribution to the effective single-particle 
potential can be written as 
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1. The thermodynamical method for transition density in 
neutron stars 
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+ 2(37r2)2/3p5/3^^^ 
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2 , (Vp,)' 
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(29) 



where the first term in the right hand side of Eq. is 
from the direct term contribution while the other terms 
arc from the exchange term contribution. As expected, 
the direct term contribution is isospin independent while 
the exchange term contribution depends on the isospin 
asymmetry as well as the density gradients. Furthermore, 
one can see that the [/-boson contribution to the single- 
particle potential depends on both the coupling constant 
g and the mass /i in a complicated way. 



D. The transition density in neutron stars 



The transition density is the baryon number density 
that separates the liquid core from the inner crust in 
neutron stars and it plays an important role in determin- 
ing the structural properties of neutron stars such as the 
crustal fraction of total moment of inertia and the mass- 
radius relations of static neutron stars. In principle, the 
transition density can be obtained from comparing rel- 
evant properties of the nonuniform solid crust and the 
uniform liquid core mainly consisting of neutrons, pro- 
tons, and electrons {npe matter). However, this is prac- 
tically very difficult since the inner crust may contain 
the so-called "nuclear pasta" with very complicated ge- 
ometries m, [tsI - IsiI . In practice, a good approximation 
is to search for the density at which the uniform liquid 
first becomes unstable against small amplitude density 
fluctuations with clusterization. This approximation has 
been shown to produce a very small error for the ac- 
tual core-crust transition density and it would yield the 
exact transition density for a second-order phase transi- 
tion [ill, l82| - [8^ . So far, several such methods including 
the thermodynamical method 14^1851 ISo l , th e dynamical 
curvature matrix method [38l - t4ll [83 . [871 [881 1 , the Vlasov 



equation metho d I89l-l93i. and the random phase approx- 
imation (RPA) [8491. 94 1 have been applied extensively 
in the literature. Here, we briefly introduce the ther- 
modynamical method, the dynamical curvature matrix 
method and the Vlasov equation method, which will be 
used to calculate the transition density in this work. 



In the thermodynamical method, the system is re- 
quired to obey the following intrinsic stability condi- 
tion [H,!!!!!! 



dP 

dv 

dp,. 



■np 



> 0, 



> 0, 



(30) 
(31) 



where the P = Pb + Pe the total pressure of the npe 
matter system with P^ and Pe denoting the contributions 
from baryons and electrons respectively, and the v and qc 
are the volume and charge per baryon number. The 
is defined as the chemical potential difference between 
neutrons and protons, i.e.. 



Pnp — Pn pp • 



(32) 

The conditions of Eq. ((30)) and Eq. ((3T|) are equivalent 
to require the convexity of the energy per particle in the 
single phase [H, [sl] by ignoring the finite size effects due 
to surface and Coulomb energies [9^. In fact, Eq. ([50)) is 
simply the well-known mechanical stability condition of 
the system at a fixed /i„p, which ensures that any local 
density fluctuation will not diverge. On the other hand, 
Eq. ([3T|) is the charge or chemical stability condition of 
the system at a flxed density. It means that any local 
charge variation violating the charge neutrality condition 
will not diverge. 

The pressure Pe is only a function of the chemical po- 
tential difference pnp by assuming the /^-equilibrium con- 
dition is satisfied, i.e., pnp = Pe- By using the relation 
^^aa'"^"^ ^ ~^^np with Eb{p, Xp) being energy per baryon 
from the baryons in the /3-equilibrium neutron star mat- 
ter and Xp = Pp/p, and treating the electrons as free 
Fermi gas, one can show [9^ that the thermodynamical 
relations Eq. ([30]) and Eq. ([31]) are actually equivalent to 
the following condition 



Vthe 



2p 



dEb{p,Xp) 29'^^b{p,Xp) 



dp 



+ p- 



d'^Eb{p,X p) \^ ^d'^Eb{p,Xp) 
dpdxp 



P] I- 



dx^ 



> 0, (33) 



which determines the thermodynamical instability region 
of the /3-equilibrium neutron star matter. The baryon 
number density that violates the condition Eq. p3|) then 
corresponds to the core-crust transition density in neu- 
tron stars for the thermodynamical method. 



2. The curvature matrix method for transition density in 
neutron stars 

In the curvature matrix method, the instability region 
of homogeneous nuclear matters against clusterization is 
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determined by introducing a finite-size spatially periodic 
density fluctuation Sp to the system and then examin- 
ing how the system free energy varies with the fluctua- 
tion fssj. The fluctuation will affect the three compo- 
nents of homogeneous nuclear matter (neutrons, protons 
and electrons) independently when assuming it occurs 
only on finite microscopic scale in the /3-equilibrium nu- 



clear matter as 

Pq=Pq+ 6pg, (34) 

with q = n, p, e. Then the free energy / at each point of 
density Pq{r) = + 5pq can be expressed as 



q,q' —n,p,e 



2 \dpqdpq' J ^ 



SpqSpq, 



(35) 



For a stable homogeneous npe matter system, the first-order term (^Sp~^ ^EH) must equal to zero and a density 

fiuctuation 6pq should lead to an increasing of the free energy, which is equivalent to require the second-order term 
in Eq. ([55]) to be positive for any density fluctuation Spq. This can be ensured by the positive definiteness of the 
following curvature matrix 



/ 



C, 



CM 



\ dpadp 
dp 





„ dpcdpp 

dpp 
dpp 
dpp 





dpndpe 

dppdpc 



9pe 




(36) 



where k is the wave vector of the spatially periodic den- 
sity fluctuations and the effective density-gradient coef- 
ficients are defined as 

3 

DnniDpp) = —[tl{l~Xi)^t2{l+X2)] 

+ T^(^^Mrt+P«(p)a^j> (37) 

Dnp{Dpn) = ^[3ti(2 + xi) -t2(2 + a;2)] 

- ^[ti(2 + xi)-H<2(2-f xa)]. (38) 

From Eq. ((36|) . one can see that the matrix C^j^j includes 
four parts. The bulk part, i.e., the first term in the right 
hand side of Eq. ([55]) is k independent but the following 
3 parts are all k dependent due to the density-gradient 
terms in Eq. (jlSp . the direct term contribution of the 
finite range interaction from the [/-boson exchange, and 
the Coulomb interaction, respectively. 

The matrix Cqj^^ is defined for each point {pn,Pp,Pe), 
and the sign of its three eigenvalues determines the sign 
of the second-order term in Eq. (|35p . namely, only if all 



the eigenvalues of the matrix are positive, the free en- 
ergy of the system will remain the minimum value and 
the nuclear system will be stable for all the density fluc- 
tuations. So the baryon number density violating the 
positive definiteness of the matrix C^j^^ corresponds to 
the core-crust transition density in neutron stars for the 
curvature matrix method. 



3. The Vlasov equation method for transition density in 
neutron stars 



To determine the core-crust transition density in neu- 
tron stars within the Vlasov equation method, we include 
here for completeness a brief description for the method 
(See, e.g., Ref. [H IH for the dctafls). For a /3-stable 
and electrically neutral npe matter, the Vlasov equation 
can be expressed as 



dfqjR.P^t) 
dt 



'Vq-Vgfq (i?, p,t)~V^Uq- Vpfq {R,p,t)=Q 

(39) 



in terms of the semi-classical Wigner function for particle 
type q = n,p,e, i.e., 

X (R+^-,t)e^P-^^d\, (40) 



with Sq = Lo/kVq and = Pq/{m*'^ + Pq'^Y^'^ being 
the Fermi velocity. The momentum integration can be 
evaluated approximately as 



^1 (-m.'* 



1 
2^2 



dea 



F * 
Pq mq 

27r2 



(46) 



where (jjqi is the wave function of z-th particle of type q, R 
and r are defined as the same as R and r in the previous 
section. In Eq. ([M)) . Vq = p/{m*^ +p2)i/2 denotes the 
particle velocity and m* — m^. 

The density fluctuation due to a collective mode with 
frequency lu and wavcvcctor k in nuclear matter can be 
studied through following the standard procedure ^ by 
writing 



fq{R,P,t)^f^ip)+6fqiR,p,t) 



with 



we obtain the Vlasov equation 



6pq « XqLq J^^^'^'^ ' 



where Lq is the usual Lindhard function 



(41) 



6fq{R,p,t)^Sfq{p)e-^-'+^'-^. (42) 
By expressing pq{R,t) = Pq + Spq{R,t) with 

^ '' Sfq{R,p,t)d^p, (43) 



(44) 



for q = n,p with « Pq /2m* and 



J4_ 

27r2' 



(47) 



for electrons, where p,e « is the electron chemical 
potential. 

For protons, there are additional direct and exchange 
Coulomb contributions to the factor j^5pqi in 



Eq. p4)) given, respectively, by 



CD 



SU, 



47re^ 



{Spp - Spe), 



1 2 



1/3 



(48) 



(49) 



For electrons, there are only direct and exchange 
Coulomb contributions to ^^Spq'. 

After linearizing the Vlasov equation, we can reexprcss 
Eq. as a function of the collective density fluctuation 



(50) 



;i cos6lrf(cos6i) ^ , fsq + l\ 
^^-/_, .,-cos/ --^ + -^^H^)' ^''^ 

I 



with 



Y T ^ 



XpLp 



dUr, 



XnLn 






XpLp 



X T ^ 
^P^P dp. 



X„L, 






XpLf 



R 

dpp 





r dUp 



dp„ 





XeL. 
au„ 

dpp 
dUp 
dpp 
dlh 
dpp 





dUp 
dp. 

















1 






^0 








(51) 



In Ec^. (I5ip . the C/^ is the single-particle potential in the npe system and the matrix 
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has the same fc-dependent terms as in Eq. ((36)) . 

The determinant |Cy^| = determines the dispersion 
relation u}{k) of the collective density fluctuation which 
also determines the non-trivial solutions of Eq. ([SUI . The 
transition density in neutron stars is the density at which 
the frequency ui becomes imaginary, leading to that the 
collective density fluctuation would grow exponentially 
and thus the instability of the neutron star matter oc- 
cur. To determine the condition for this to occur, we 
let Sq = —ivq {vq > 0) and rewrite the Lindhard func- 
tion as Lq = —2 + 2;^^ arctan(l/i'q). Since the values of 
Lq are in the range of —2 < Lq < 0, the critical values 
Ln = Lp = Le = —2, corresponding to Vq = 0, then de- 
termine the spinodal boundary of the system when they 
are substituted into |Cy^| = 0. The baryon number 
density that makes |Cyg| vanish then corresponds to the 
spinodal boundary in the neutron star matter or the core- 
crust transition density of neutron stars for the Vlasov 
equation method. 



III. RESULTS 

In the present work, for the Skyrme effective nucleon- 
nucleon interaction, we u se the modified Skyrme-like 
(MSL) parameter [11 [ToOl for which the 9 Skyrme inter- 
action parameters a, tQ—t^, xq — xs are obtained analyti- 
cally in terms of 9 macroscopic quantities po? ^o(Po); the 
incompressibility Kq, the isoscalar effective mass m* q, 
the isovector effective mass toJq, i?sym(po)i L, the gra- 
dient coefficient Gs, and the symmetry-gradient coeffi- 
cient Gv- In particular, the MSLO parameter set [93] is 
obtained by using the following empirical values for the 
9 macroscopic quantities: po = 0.16 fm~'^, i?o(po) = ^16 
MeV, Kq = 230 MeV, ?«* o = 0.8m, ml_Q = 0.7m, 
Ssym(po) = 30 MeV, L = 60 MeV, = 5 MeV-fm^, 
and Gs = 132 MeV-fm^. And the spin-orbit coupling 
constant Wq = 133.3 MeV •fm'"' is used to fit the neutron 
Pi/2 " P3/2 splitting in ^^O. It has been shown [93] that 
the MSLO interaction can give a good description of the 
binding energies and charge rms radii for a number of 
closed-shell or semi-closed-shell nuclei. 

For the J7-boson-nucleon coupling constant g and the 
[/-boson mass /i, their values are largely uncertain. As 
argued by Krivoruchenko efc al. (soj . in order to ensure 
that the [/-boson effects on finite nuclei should be neg- 
ligible, the Compton wavelength of the [/-boson is usu- 
ally assumed to be greater than the radius of heavy nu- 
clei, i.e., about 7 fm, leading to ^ < 30 MeV. At the 
same time, the .g^//i^ value of the [/-boson should be 
less than about 200 GeV~^, which roughly corresponds 
to the value of the ordinary vector a; boson. Otherwise, 
the [/-boson is neither weakly coupled nor light. On the 
other hand, there also exist some constraints on prop- 
erties of the [/-boson from cosmology and astrophysical 
observations. For example, the [/-boson mass fi is re- 
quired to exceed the mass of light cold dark matter, i.e.. 



MeV, to explain the excess fiux of 511 keV photons 
coming from the central region of our Galaxy observed 
by the SPI/INTEGRAL satellite [lOl|- Based on above 
discussions, in the present work we assume the [/-boson 
mass is in the range of 2 MeV < M < 30 MeV (the cor- 
responding Compton wavelength of the [/-boson is thus 
between about 7 fm and 100 fm) and the g^/fJ.^ value of 
the [/-boson satisfies g^/fi^ < 150 GeV~^. The latter is 
further consistent with existing constraints from neutron- 
proton and neutron-lead scatterings, the spectroscopy of 
antiproton atoms as well as the recently discovered new 
holder of neutron star maximum mass of 1.97 ± 0.04 Mq 
from PSR J1614-2230 [H, [11, [H, [H . 

A. Nuclear matter symmetry energy from the 
[/-boson 

As have been seen in the previous section, the direct 
term of the [/-boson contributes to the nuclear matter 
EOS only through the combination g^ / fi^, and in partic- 
ular, its contribution to the energy per nuc leon is given 
by ^g'^ / y?p. As it was emphasized by Fujii jl02j . for the 
direct term contribution, though both the coupling con- 
stant g and the mass p are small for the light and weakly 
coupled bosons, the value of the ratio / IJ? can be large. 
Therefore, the light and weakly coupled bosons can sig- 
nificantly affect the nuclear matter EOS and thus the 
properties of neutron stars (30l - [33j . On the other hand, 
the exchange term contribution to the nuclear matter 
EOS depends on both the coupling constant g and the 
mass /i in a complicated way. Furthermore, it is interest- 
ing to see that the isoscalar [/-boson will contribute to 
the nuclear matter symmetry energy due to the exchange 
term contribution though the direct term does not have 
such contribution. 

To see quantitatively how the exchange term of the 
[/-boson affects the nuclear matter EOS, we show in 
Fig. [T]the exchange term contribution of the [/-boson to 
the EOS of symmetric nuclear matter £^^ub(p) ((^) ^^"^ 
(b)) and the nuclear matter symmetry energy i?sym,UB(p) 
((c) and (d)) as functions of density with several typi- 
cal values of p. and g^ j [i^ ■ Here, the £'sym,UB(p) is ex- 
tracted from the parabolic approximation iS'sym.UBCp) ~ 
-Eub(p, (5 = 1) - £;ub(p,'5 = 0) with £'ub(p,'5) being the 
energy per nucleon from the U -boson contribution. One 
can see from Fig. [T]that both E^_^^{p) and i?sym,UB(p) 
increase with increasing values of both p and g^/p^. 
As expected, however, both the contributions of the ex- 
change term to energy per nucleon of symmetric nuclear 
matter and the symmetry energy are quite small and can 
be safely neglected compared with the direct term con- 
tribution for the values of p and g^ j p^ considered here. 
For example, even at very high baryon density such as 
p ^ 1.0 fm-3 with p = 30.0 MeV and g'^/p'^ = 150 
GeV~^ (See the right panels of Fig. [Ij, the magnitude 
of the Eqtj^{p) is only about 1.1 MeV and the magni- 
tude of the i?sym,UB(p) IS Icss than about 0.32 MeV while 



9 



w 

I 



> 



w 

I 



1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 

0.3 

0.2 

■0.1 



:(a) n= 1.97 MeV , - 

,' x200 

✓ 

- ✓ 

_/ ^ 1<200 


:(b) n = 30.0 MeV' 

✓ 


L - il^ = 150 GeV"'^ 
f gV^VSGeV"' 


/ 

' ^ " 

' y — 


1 1 1 1 1 1 1 

(c)n= 1.97 MeV 

^ - ' x200 


-(d) n = 30.0 MeV 

✓ 

/ -- — ' 







0.0 

0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 1.0 

P (fm ') 

FIG. 1: (Color online) The exchange term contribution of the 
[/-boson to the EOS of symmetric nuclear matter ((a) and 
(b)) and the symmetry energy ((c) and (d)) as functions of 
density with several typical values of [i and / . Note: The 
results with ^ = 1.97 MeV ((a) and (c)) have been rescaled 
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FIG. 2: (Color online) Transition density pt (a) and pressure 
Pt (b) in neutron stars as functions of the L parameter with 
the MSL interaction using the thermodynamical method, the 
curvature matrix method and the Vlasov equation method. 



the direct term contribution to the energy per nucleon of 
symmetric nuclear matter reaches about 576 MeV. The 
magnitudes of Eq^^{p) and -Bsym.UBCp) will further de- 
crease if smaller values of ^ and j [P' are used (See, e.g., 
the left panels of Fig. [T|). These results verify the valid- 
ity of neglecting the exchange term contribution of the 
?7-boson to nuclear matter EOS in the literature (30l - l33| . 

B. The core-crust transition density and pressure 
in neutron stars with the L*^-boson 

We now turn to the numerical results on the core-crust 
transition density and pressure in neutron stars with 
the thermodynamical approach, the curvature-matrix ap- 
proach, and the Vlasov equation approach. We note that 
both the matrices C^^.^^ in Eq. (p6)) and Cy^ in Eq. ((5T|) 
are fc-dependent, and for the curvature-matrix approach 
and the Vlasov equation approach, the core-crust tran- 
sition density corresponds to the critical baryon number 
density above which the neutron star matter is always 
stable for all possible values of k while below which one 
can always find a fc value to violate the stability condi- 
tions of the neutron star matter. On the other hand, for 
the thermodynamical method, the transition density can 
be directly obtained by solving the equation Vther = 
(See Eq. ([55)1 for the expression of Vther)- 

Theoretically it has been established that there ex- 
ists a strong correlation between the transition density 
Pt and the nuclear symmetry energy. In particular, a 
strong linear correlation between the transition density 
Pt and the slope parameter L of the nuclear symmetry 



energy has been observed in many different theoretical 
calculations [96l - [99j . To see the symmetry energy depen- 
dence of the inner edge of neutron star crusts, we show in 
Fig. [5] the transition density pt and pressure Pt in neu- 
tron stars as functions of the L parameter with MSLO 
interaction by varying individually L using the thermo- 
dynamical method, the curvature matrix method and the 
Vlasov equation method. When varying individually the 
L parameter, we keep all other macroscopic quantities 
Po, Eo{po)., Kq, m*_o, TO* 0, £^sym(po), Gs, Gy , and Wo 
at their default values in MSLO. It should be noted that 
the original agreement of MSLO with the experimental 
data of binding energies or charge radii of finite nuclei 
essentially still holds with the individual change of the L 
parameter. 

For the results shown in Fig. [21 the [/-boson contri- 
butions are not considered. It is seen that all the three 
methods give similar results for the L dependence of the 
transition density pt and pressure Pt with the curvature- 
matrix method giving slightly smaller values of pt and Pt 
than the thermodynamical method while slightly higher 
values of pt and Pt than the Vlasov equation method for 
a fixed value of L. The nonmonotonous variation of the L 
dependence of the transition pressure Pt in the thermody- 
namical method is due to the fact that the Pt is a compli- 
cated functions of pt, L, and the isospin asymmetry 6t at 
the transition density (See, e.g., [l^). The smaller values 
of Pt from the curvature-matrix method than from the 
thermodynamical method implies that the density gradi- 
ent terms and Coulomb term considered in the former can 
make the neutron star matter more stable, which are con- 
sistent with the results in Ref. [9^ (The curvature- 
matrix method is called dynamical method there). On 
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FIG. 3: (Color online) The / dependence of the tran- 
sition density pt and pressure Pt in neutron stars from the 
thermodynamical method, the curvature matrix method and 
the Vlasov equation method with the MSLO interaction for 
= 1.97 MeV and 30.0 MeV, respectively. 



the other hand, the slightly smaller values predicted by 
the Vlasov equation method than the curvature-matrix 
method is due to the quantum effects considered in the 
Vlasov equation method, indicating that the quantum 
effects will make the neutron star matter more stable as 
expected. 

To see how the light and weakly coupled vector U- 
boson affects the inner edge of neutron star crusts, we 
show in Fig. [3] the / IJ? dependence of the transi- 
tion density pt and pressure Pt in neutron stars from 
the MSLO interaction by including the J7-boson with 
/i = 1.97 MeV and 30.0 MeV, respectively, using the ther- 
modynamical method, the curvature matrix method and 
the Vlasov equation method. One can see clearly that 
while the curvature matrix method and the Vlasov equa- 
tion method predict very similar results for the 
dependence of pt and Pt, the thermodynamical method 
predicts very different results with pt and Pt decreasing 
very quickly as the j p^ increases. In particular, we 
find that for j p^ > 20 GeV~^, the neutron star matter 
is always stable and the transition density does not ex- 
ist in the thermodynamical method. This is due to the 
fact that the fc-dependent terms in Eq. ((36)) and Eq. ([ST]) 
originating from the finite-range interaction and density- 
gradient contributions play an important role in deter- 
mining the transition density when the ?7-boson is con- 
sidered. Therefore, the thermodynamical method which 
ignores the fc-dependent terms would be no longer appro- 
priate to determine the inner edge of neutron star crusts 
when the f/-boson is taken into account. 

It is interesting to see from Fig. [3] that for the cur- 
vature matrix method and the Vlasov equation method, 
the effects of the U -boson on the transition density pt and 
pressure Pt depend on not only the ratio g^ / p^ but also 
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FIG. 4: (Color online) The l/jj, dependence of the transition 
density pt and pressure Pt in neutron stars from the the cur- 
vature matrix method and the Vlasov equation method with 
the MSLO interaction for g'^/p'^ = 75 GeV"^. 



the [/-boson mass p. In particular, for a heavier [/-boson 
(e.g., p = 30.0 MeV), the pt decreases significantly with 
increment of g'^/p^ from to 150 GeV~^. On the other 
hand, for a very light [/-boson (e.g., p = 1.97 MeV), 
the transition density pt exhibits very weak dependence 
on g^/p^- These features imply that for a fixed value 
oi g^ / p^, a heavier [/-boson can make the neutron star 
matter more stable while a very light [/-boson essentially 
has no influence on the transition density. 

As shown in Fig. [3l for the curvature matrix method 
and the Vlasov equation method, although the transition 
density pt displays a somewhat complicated relationship 
with the properties of the [/-boson, the transition pres- 
sure Pt simply increases with the ratio g^ / p^ whether the 
[/-boson mass is light or heavy. Furthermore, it is inter- 
esting to see from Fig. [3] (b) that the transition pressure 
Pt increases almost linearly with g^ / p^ for a very light 
[/-boson (i.e., p = 1.97 MeV), while it exhibits much 
slower increment with g'^/p'^ for a heavier [/-boson (e.g., 
p = 30.0 MeV) as shown in Fig.[3](d). The linear correla- 
tion between Pt and g^/p^ for p = 1.97 MeV observed in 
Fig. [3] (b) is easily understood since the [/-boson contri- 
bution to the pressure Pt, ub is dominated by the direct 
term contribution, i.e., Pt, ub ~ ^9^/f^^Pt' if M is inde- 
pendent of the density as we assume in the present work, 
and the pt remains approximately a constant value for 
p = 1.97 MeV when g'^/p'^ varies from to 150 GeV"^ 
as shown Fig. [3] (a). In the case of ^ = 30.0 MeV, the 
Pt decreases significantly with the increment of g^/p^ as 
shown Fig. [3] (c), which leads to the transition pressure 
Pt displays much slower increment with g^ j p? as shown 
in Fig. 13] (d). 
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In order to see the [/-boson mass dependence of pt and 
Pt at a fixed value of j [j? , we display in Fig. 0] the 1/^ 
dependence of the transition density pt and pressure Pt 
in neutron stars from the MSLO interaction by including 
the [/-boson contribution with j = 75 GeV~^ using 
the curvature matrix method and the Vlasov equation 
method. It is seen that the two methods predict very sim- 
ilar 1/ p dependence of the transition density pt and pres- 
sure Pt with the Vlasov equation method giving smaller 
values. It is interesting to see that pt becomes sensitive 
to p when the [/-boson mass p is larger than about 2 
MeV although the [/-boson almost has no influence on 
the transition density pt if its mass p is less than about 2 
MeV. The transition pressure Pt displays similar 1/p de- 
pendence as the Pt due to the relation Pt,uB ~ ^g^/p^Pt- 
Therefore, these results demonstrate that the transition 
density pt in neutron stars can be sensitive to the value 
of both p and g^ j p? ^ and any experimental or observa- 
tional constraints on pt may put important limits on p 
and g^ I p} , or equivalently on p and g. 

As have been shown above, the exchange term contri- 
bution of the [/-boson to the nuclear matter EOS can be 
safely neglected. It is thus interesting to see how the ex- 
change term affects the inner edge of neutron star crusts. 
To check this point, we show in Fig. [S] the Xjp depen- 
dence of the transition density pt in neutron stars with 
the MSLO interaction for g^ j p^ = 75 GeV"^ from the 
Vlasov equation method together with that by neglecting 
the direct term contribution or neglecting the exchange 
term contribution. It is seen that the exchange term 
contribution of the [/-boson has very small influence on 
the transition density pt, especially for the light mass [/- 
boson. On the other hand, the direct term signiflcantly 



affects the Xjp dependence of the transition density pt. 
Particularly, neglecting the direct term contribution of 
the [/-boson leads to very weak Xjp dependence of the 
transition density pt, implying that the observed strong 
Xj p dependence of the transition density pt is essentially 
due to the direct term contribution which produces a fc- 
depcndent term in the matrix (|36p or (|5ip . i.e., the third 
term in the right hand side of Eq. (pS)) . We note that 
using the curvature matrix method leads to the same 
conclusion. 



C. The mass-radius relation and crustal fraction of 
moment of inertia for static neutron stars with the 
[/-boson 

As we have showed above, the [/-boson may have sig- 
nificant influence on the nuclear matter EOS and the 
inner edge of neutron star crusts. Here we investigate 
effects of the [/-boson on the global properties of static 
neutron stars. To calculate the global properties, such 
as the mass-radius relation and crustal fraction of mo- 
ment of inertia, of static neutron stars, one needs the 
EOS of neutron star matter over a broad density region 
ranging from the center to the surface of neutron stars. 
Besides the possible appearance of nuclear pasta in the 
inner crust, various phase transitions and non-nucleonic 
degrees of freedom may appear in the core of neutron 
stars. In this work, we restrict ourselves to the simplest 
and traditional model, and make the minimum assump- 
tion that the core of neutron stars contains the uniform 
/^-stable and electrically neutral npep matter only and 
there is no phase transition. 

Generally, a typical neutron star contains the liquid 
core, inner crust and outer crust from the center to sur- 
face. For the liquid core we use the EOS of wpep mat- 
ter from SHF calculations including the [/-boson con- 
tributions to the nuclear EOS. For the Skyrme effective 
nucleon-nucleon interaction, the MSL interaction with a 
soft symmetry energy of L = 30 MeV is used. We note 
that the MSL interaction with L = 30 MeV predicts a 
n-pep matter EOS very similar to the more sophisticated 
EOS containing nucleons, hyperons and quark degrees of 
freedom [33|. In the present work, the [/-boson contri- 
bution to the EOS of liquid core includes both the di- 
rect term contribution (i.e., Eq. (HH) and the exchange 
term contribution (i.e., Eq. (j23p without the density gra- 
dient terms) although the latter is negligible. In par- 
ticular, the fractions of neutrons, protons, electrons and 
muons in the neutron star matter are obtained from self- 
consistently solving the set of equations for /3-stable con- 
dition (i.e., pnp = Me = P^i) and charge neutral con- 
dition(i.e., Pp = Pe + Pp) by considering the [/-boson 
exchange term contribution to the chemical potential of 
neutrons and protons. It should be noted that the U- 
boson direct term does not change the neutron and pro- 
ton chemical potential difference pnp as it contributes 
equally to the single-particle potential of neutrons and 
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protons, and thus the chemical compositions of the neu- 
tron star matter will not change if only the J7-boson di- 
rect term contribution is considered as pointed out in 
previous work (30l - [33l |. In this way, the contributions of 
U -boson to the energy density, the pressure, the nucleon 
effective masses and chemical potentials, are considered 
self-consistently in the neutron star matter calculations 
for the liquid core. 

In the inner crust with densities between pout and pt 
where the nuclear pastas may exist, because of our poor 
knowledge about its EOS from first principle, following 
Carriere et al. [s^l (See also Ref. [9g) we construct its 
EOS according to 



(53) 



This polytropic form with an index of 4/3 has been found 
to be a good approximation to the crust EOS [i^. |45|. 
The Pout = 2.46 x 10^* fm^'^ is the density separating 
the inner from the outer crust. The constants a and b 
are then determined by 



P c 

^ out^t 



4/3 



p 4/3 



g4/3 



Pt - Po . 



4/3 4/3 ' 



(54) 
(55) 



where Pt, et and Pout, £out are the pressure and energy 
density at pt and Pout, respectively. In the outer crust 
with 6.93 X 10-^3 fm-^ < p < pout, we use the EOS of 



103|, and in the region of 4.73 x IQ-^^ fm"-^ < 



BPS . , 

p <6.93 X lb~^^ fm""* we use the EOS of Feynman- 
Metropolis- Teller (FMT) For the t/-boson contri- 

bution to the EOS of neutron star crusts and surface, we 
add the energy density and pressure from the J7-boson 
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FIG. 7: (Color online) The mass-radius relation of static neu- 
tron stars using the MSL interaction with L = 30 MeV and 
different values of / for ^ = 1.97 MeV (solid lines) and 
30.0 MeV (dashed lines), respectively. The results with pt 
from both the Vlasov equation method (thick lines) and the 
curvature matrix method (thin lines) are included for com- 
parison. The shaded band represents the latest new holder of 
the maximum mass of neutron stars of 1.97 ± 0.04Mq from 
PSR J1614-2230 [H. 



direct term contribution, i.e., euB = Pub ~ ^S^/a'^P^ to 
the corresponding parts since the exchange term contri- 
bution is negligible as shown in the above. 

As an example, we show in Fig.|6]the EOS for different 
parts of a neutron star. As we have discussed earlier, the 
transition density pt is obtained by studying the onset 
of instabilities in the liquid core, namely it is the critical 
density below which small density fiuctuations will grow 
exponentially. Therefore, the transition density pt and 
the EOS of the liquid core are obtained self-consistently 
from the same interaction and in this sense they are on 
the same footing. We use in Fig.[n]the pt obtained within 
the Vlasov equation method using the full EOS with the 
MSL interaction of L = 30 MeV for different values of 
g'^/p? with p. = 1.97 MeV. Using the above EOS for 
different parts of the neutron star, the radial distribution 
of the total energy density and the pressure in neutron 
stars is continuous, but the derivative of the pressure is 
not continuous at pt and pout- It is seen that the EOS's 
of the inner crust and the liquid core are quite different 
for different values of j pi . Interestingly, one can see 
that because the transition density pt displays very weak 
dependence on g^ j pi for the very light J7-boson mass 
(1.97 MeV here), the pt (and thus the corresponding £t) 
has essentially the same value for different g'^ / pi, and 
the difference of the EOS for different values of g'^ / pi 
observed in Fig. [HI is essentially due to the variation of Pt 
with g^ / pi because of the relation Pt^uB ~ \g'^ / pipl- 

Using the EOS constructed above, one can solve the 
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Tolman-Oppenheimer-Volkoff (TOV) equations to obtain 
the mass-radius relations and the results are shown in 
Fig. [71 Indicated by the shaded band in Fig. [7] is the 
latest new holder of the maximum mass of neutron stars 
of 1.97 ± 0.04Mq from PSR J1614-2230 H^. For MSL 
interaction with a soft symmetry energy (L = 30 MeV) 
without considering the t/-boson contribution, the neu- 
tron star mass M decreases quickly with increasing ra- 
dius R and the maximum mass is about 1.47 Mq, which 
is significantly less than the observed maximum neutron 
star mass of 1.97±O.O4M0. On the other hand, the neu- 
tron star mass can be enhanced strongly if the effects of 
[/-boson are considered. In particular, a larger value of 
g'^ / /i^ leads to a larger neutron star mass at a fixed radius 
since the nuclear EOS is increasingly stiffened with incre- 
ment of I [i^ as shown in Fig. [5] Especially, the neutron 
star maximum mass can reach IXsIMq with j [i^ = 75 
GeV~^, and the corresponding radius of the maximum 
mass neutron star is about 14.8 km. 

To see more clearly the [/-boson effects on the prop- 
erties of neutron stars, we include in Fig. [7] the mass- 
radius relations for different values of g^ j [i^ with two 
values of the [/-boson mass, i.e., /i = 1.97 and 30.0 MeV, 
respectively. Furthermore, the results with pt from both 
the Vlasov equation method and the curvature matrix 
method are included for comparison although the former 
is believed to be more realistic. As expected, for the 
larger value of the [/-boson mass, a smaller value of the 
neutron star mass M at a fixed radius R is obtained due 
to the smaller transition density and pressure obtained 
for a larger p as shown in the previous section. This 
[/-boson mass effect on the neuron star mass-radius re- 
lation will become more pronounced for a larger value of 
g^ / p^ (e.g., g^ I [i^ = 75 GeV^^). In addition, one can see 
that using the pt from the curvature matrix method pre- 
dicts a little larger neutron star mass M at a fixed radius 
R than using the pt from the Vlasov equation method. 
However, the maximum mass of the neutron stars is al- 
most the same for different values of \x using the pt from 
either the curvature matrix method or the Vlasov equa- 
tion method. These results indicate that, essentially re- 
gardless of the values of the [/-boson mass and the pt, a 
value of g"^ / p? = 75 GeV~^ can reasonably describe the 
latest new holder of the neutron star maximum mass of 
1.97 ± 0.04AfQ from PSR J1614-2230 H^. 

The crustal fraction of total moment of inertia of a 
static neutron star, A///, is a particularly interesting 
quantity as it can be inferred from observations of pulsar 
glitches, i.e., the occasional disruptions of the otherwise 
extremely regular pulsations from magnetized, rotating 
neutron stars. Furthermore, as was stressed by Lattimer 
and Prakash [i^l, the AI/I depends sensitively on the pt 
and Pt which are essentially determined by the EOS of 
asymmetric nuclear matter at subsaturation densities as 
shown above, but there is no explicit dependence upon 
the EOS of neutron star matter at higher densities. These 
features imply that the AI/I may provide a good probe 
for properties of the [/-boson. To illustrate the [/-boson 
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FIG. 8: (Color online) The crustal fraction of total moment of 
inertia AI/I as a function of the neutron star mass for static 
neutron stars using the MSL interaction with L = 30 MeV for 
g'/p^ = 75 GeV"^ with p = 1.97 MeV (solid lines) and 30.0 
MeV (dashed lines), respectively. The results with pt from 
both the Vlasov equation method (thick lines) and the curva- 
ture matrix method (thin lines) are included for comparison. 
The lower limit AJ/J = 0.014 of the observation constraint 
for the Vela pulsar [i^] is also indicated. 



mass dependence of the AI/I, we show in Fig. [8] the 
AI/I as a function of the neutron star mass for static 
neutron stars using the MSL interaction with L = 30 
MeV for g^/p^ = 75 GeV"^ with p = 1.97 and 30.0 MeV, 
respectively. Furthermore, the results with pt obtained 
from both the Vlasov equation method and the curvature 
matrix method are included for comparison. The AI/I 
is obtained here from direct numerical calculations as in 
Ref. [9^. As expected, one can see that the AI/I indeed 
exhibits a clear sensitivity to the [/-boson mass with a 
heavier [/-boson mass giving a smaller AI/I for a fixed 
neutron star mass. In addition, one can see that using 
the Pt obtained from the curvature matrix method gives 
a little larger value for AI/I at a fixed neutron star mass 
M than using the pt from the Vlasov equation method. 
Empirically, the crustal fraction of total moment of iner- 
tia has been constrained as AI/I > 0.014 from studying 
the glitches of the Vela pulsar [l^ and the lower limit 
AI/I = 0.014 of the constraint is also indicated in Fig. 
m It is seen from Fig. [8] that the calculated results of 
AI/I with the two values of = 1.97 and 30.0 MeV 
using the pt obtained from either the curvature matrix 
method or the Vlasov equation method are all consistent 
with the observation constraint of AI/I > 0.014 for the 
Vela pulsar. 

Based on the results above, we conclude that the vec- 
tor [/-boson can significantly stiffen the nuclear matter 
EOS and thus enhance strongly the (maximum) mass of 
neutron stars. Furthermore, as the pt and Pt are sensi- 
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tive to both j [i^ and [i if the f/-boson mass is larger 
than about 2 MeV, both / jJ? and /z can thus affect the 
mass-radius relation of neutron stars although the max- 
imum mass of neutron stars is essentially independent 
of the J7-boson mass In addition, our results demon- 
strate that the crustal fraction of total moment of inertia 
A/// may depend sensitively on the J7-boson mass /i for 
a fixed value of 

IV. SUMMARY 

Using the thermodynamical approach, the curvature 
matrix approach and the Vlasov equation approach with 
the Skyrme effective nucleon-nucleon interaction, we have 
investigated effects of the light vector gauge t/-boson, 
that is weakly coupled to nucleons, on the core-crust 
transition density pt and pressure Pt of neutron stars. For 
the exchange term contribution of the ?7-boson, we have 
applied the density matrix expansion approach, which 
automatically leads to the density gradient terms in the 
single nucleon potential and nuclear energy density func- 
tional. Our results have demonstrated that the exchange 
term contribution to energy per nucleon of symmetric nu- 
clear matter and the symmetry energy is quite small and 
can be safely neglected compared with the direct term 
contribution for the parameter range of ^ and / p^ con- 
sidered in this work, verifying the validity of neglecting 
the [/-boson exchange term contribution to the nuclear 
matter EOS as have been done in the literature. Fur- 
thermore, the exchange term has also been found to have 
negligible influence on the core-crust transition density pt 



and pressure Pt, especially for very light /7-boson. 

Interestingly, our results have shown that the pt and Pt 
depend on not only the ratio of coupling strength to mass 
squared of the [/-boson g^ j p^ but also its mass p. The 
[/ -boson mass dependence of pt and Pt is due to the finite 
range interaction from the [/-boson exchange. Especially, 
we have found that the pt and Pt will be sensitive to both 

I p^ and p if the [/-boson mass p is larger than about 2 
MeV, and both g^ j p^ and p can have significant influence 
on the mass-radius relation and the crustal fraction of 
total moment of inertia of neutron stars. Therefore, our 
results presented in this work have demonstrated that as- 
trophysical observations on neutron star structures, such 
as the mass-radius relation and the crustal fraction of 
total moment of inertia from pulsar glitches, can be po- 
tentially useful to constrain properties of the [/-boson, 
e.g., its mass p and the coupling constant g to nucleons. 



ACKNOWLEDGMENTS 

The authors would like to thank W.Z. Jiang, B.A. Li, 
D.H. Wen, and J. Xu for useful discussions. This work 
was supported in part by the National Natural Science 
Foundation of China under Grant Nos. 10975097 and 
11135011, Shanghai Rising-Star Program under grant 
No. 11QH14G11GG, "Shu Guang" project supported by 
Shanghai Municipal Education Commission and Shang- 
hai Education Development Foundation, and the Na- 
tional Basic Research Program of China (973 Program) 
under Contract No. 2007CB815004. 



[1] P. Fayet, Phys. Lett. 95B, 285 (1980); Nucl. Phys. B 

187, 184 (1981). 
[2] P. Jean et al, A&A 407, L55 (2003). 
[3] C. Boehm, D. Hooper, J. Silk, M. Casse, and J. Paul, 

Phys. Rev. Lett. 92, 101301 (2004). 
[4] C. Boehm, P. Fayet, and J. Silk, Phys. Rev. D 69, 

101302(R) (2004). 
[5] C. Boehm and P. Fayet, Nucl. Phys. B 683, 291 (2004). 
[6] Y. Fujii, Nature (London), Phys. Sci. 234, 5 (1971). 
[7] N. Arkani-Hamed et al., Phys Lett. B429, 263 (1998); 

Phys. Rev. D59, 086004 (1999). 
[8] E. Fischbach and C.L. Talmadge, The Search for Non- 
Newtonian Gravity, Springer- Verlag, New York, Inc. 
(1999), ISBN 0-387-98490-9. 
[9] R. Pease, Nature 411, 986 (2001). 
[10] CD. Hoyle, Nature 421, 899 (2003). 
[11] J.C. Long et al., Nature 421, 922 (2003). 
[12] E. G. Adelberger et al, Annu. Rev. Nucl. Part. Sci. 53, 

77 (2003); Prog. Part. Nucl. Phys. 62, 102 (2009). 
[13] J. P. Uzan, Rev. Mod. Phys. 75, 403 (2003). 
[14] R.S. Decca et al., Phys. Rev. Lett. 94, 240401 (2005). 
[15] S. Reynaud et al. Int. J. Mod. Phys. A20, 2294 (2005). 
[16] Y.N. Pokotilovski, Phys. At. Nucl. 69, 924 (2006). 
[17] D.J. Kapner et al, Pliys. Rev. Lett. 98, 021101 (2007). 



[18] V.V. Nesvizhevsky et al, Phys. Rev. D 77, 034020 
(2008). 

[19] Y. Kamyshkov, J. Tithof, and M. Vysotsky, Phys. Rev. 

D 78, 114029 (2008). 
[20] M. Azam, M. Sami, C.S. Unnikrishnan, and T. Shi- 

romizu, Phys. Rev. D 77, 101101 (2008). 
[21] R.D. Newman, E.G. Berg, and P.E. Boynton, Space Sci. 

Rev. 148, 175 (2009). 
[22] A. A. Geraci et al, Phys. Rev. Lett. 105, 101101 (2010). 
[23] D. M. Lucchesi and R. Person, Phys. Rev. Lett. 105, 

231103 (2010). 
[24] P. Fayet, Phys. Rev. D 75, 115017 (2007). 
[25] S.H. Zhou, Phys. Rev. D 75, 115004 (2007). 
[26] C.H. Chen, G.Q. Geng, and G.W. Kao, Phys. Lett. 

B663, 400 (2008). 
[27] P. Fayet, Phys. Lett. B675, 267 (2009). 
[28] R. Barbieri and T.E.O. Ericson, Phys. Lett. 57B, 270 

(1975). 

[29] V.V. Nesvizhevsky and K.V. Protasov, Glass. Quan. 

Grav. 21, 4557 (2004). 
[30] M.I. Krivoruchenko, F. Simkovic, and A. Faessler, Phys. 

Rev. D 79, 125023 (2009). 
[31] D.H. Wen, B.A. Li, and L.W. Ghen, Phys. Rev. Lett. 

103, 211102 (2009). 



15 



[32] D.R. Zhang, P.L. Yin, W. Wang, Q.C. Wang, and W.Z. [ 

Jiang, Phys. Rev. C 83, 035801 (2011)^ 

D.H. Wen, B.A. Li, and L.W. Chen, [arXi^: 1101. 15041 [69 
W. Reisdorf et al. (FOPI Collaboration), Nucl. Phys. 
A781, 459 (2007). 

Z.G. Xiao, B.A. Li, L.W. Chen, G.C. Yong, and M. [70 
Zhang, Phys. Rev. Lett. 102, 062502 (2009). 
P.B. Demorest, T. Pennucci, S.M. Ransom, M.S.E. 
Roberts, and J.W.T. Hessels, Nature 467, 1081 (2010). [71 
J. Xu, L.W. Chen, CM. Ko, and B.A. Li, Phys. Rev. C 
81, 055803 (2010). 

G. Baym, C. Pethick, and P. Sutherland, Astrophys. J. 
170, 299 (1971). [72 
G. Baym, H.A. Bethe, and C.J. Pethick, Nucl. Phys. 
A175^ 225 (1971). [73; 
C.J. Pethick and D.G. RavenhaU, Ann. Rev. Nucl. Part. [74 
Sci. 45, 429 (1995). [75; 
C.J. Pethick, D.G. RavenhaU, and CP. Lorenz, Nucl. [76 
Phys. A584, 675 (1995). 

J.M. Lattimer and M. Prakash, Phys. Rep. 333-334, [77 
121 (2000); Astrophys. J. 550, 426 (2001). 

J.M. Lattimer and M. Prakash, Phys. Rep. 442, 109 [78 

(2007) . 

A. W. Steiner, M. Prakash, J.M. Lattimer, and P.J. Ellis, [79' 
Phys. Rep. 410, 325 (2005). [80 

B. Link, R.I. Epstein, and J.M. Lattimer, Phys. Rev. [81 
Lett. 83, 3362 (1999). [82 

C. J. Horowitz et al, Phys. Rev. C 69, 045804 (2004); 

C. J. Horowitz et al., Phys. Rev. C 70, 065806 (2004) [83; 

A. Burrows, S. Reddy, and T.A. Thompson, Nucl. Phys. [84 
A777, 356 (2006). 

B. J. Owen, Phys. Rev. Lett. 95, 211101 (2005). [85; 
S.B. Ruster, M. Hempel, and J. Schaffner-Bielich, Phys. 
Rev. C 73, 035804 (2006). [ 
J.W. Negele and D. Vautherin, Phys. Rev. C 5, 1472 
(1972); 11, 1031 (1975). [87 
J. Xu and CM. Ko, Phys. Rev. C 82, 044311 (2010). 

B. A. Li, L.W. Chen, and CM. Ko, Phys. Rep. 464, 113 [ 

(2008) . 

D. H. Youngblood, H.L. Clark, and Y.-W. Lui, Phys. [89 
Rev. Lett. 82, 691 (1999). 

Y.-W. Lui, D.H. Youngblood, Y. Tokimoto, H.L. Clark, [90 
and B. John, Phys. Rev. C 70, 014307 (2004). 
Z.Y. Ma et al, Nucl. Phys. A703, 222 (2002). [91 
D. Vretenar, T. Niksic, and P. Ring, Phys. Rev. C 68, 
024310 (2003). [92; 
G. Colo, N. Van Giai, J. Meyer, K. Bennaceur, and P. 
Bonche, Phys. Rev. C 70, 024307 (2004). [93; 
S. Shlomo, V.M. Kolomietz, and G Colo, Eur. Phys. J. 
A 30, 23 (2006). ' [94 

T. Li et al, Phys. Rev. Lett. 99, 162503 (2007). 
U. Garg et al., Nucl. Phys. A78 8, 36 (2007). 
G. Colo, 2009. larXiv:0902.3739l ^l [nucl-th]. [95 
L.W. Chen and J .Z. Gu, J. Phys. G 39, 035104 (2012) [96 
|arXiv:1104.5407] . 

P. Danielewicz, R. Lacey, and W.G. Lynch, Science 298, [97 
1592 (2002). 

J. Aichelin and CM. Ko, Phys. Rev. Lett. 55, 2661 [98 
(1985). 

C. Fuchs, Prog. Part. Nucl. Phys. 56, 1 (2006). [99; 
W.D. Myers and W.J. Swiatecki, Nucl. Phys. 81, 1 
(1966). [100 

[67] K. Pomorski and J. Dudek, Phys. Rev. C 67, 044316 

(2003). [101 



V. Baran, M. Colonna, V. Greco, and M. Di Toro, Phys. 
Rep. 410, 335 (2005). 

L.W. Chen, CM. Ko, and B.A. Li, Phys. Rev. Lett. 94, 
032701 (2005); Phys. Rev. C 72, 064309 (2005); B.A. 
Li and L.W. Chen, Phys. Rev. C 72, 064611 (2005). 
M.B. Tsang, Y. Zhang, P. Danielewicz, M. Famiano, Z. 
Li, W. G. Lynch, and A. W. Steiner, Phys. Rev. Lett. 
102, 122701 (2009). 

M. Centelles, X. Roca-Maza, X. Viiias, and M. Warda, 
Phys. Rev. Lett 102, 122502 (2009); M. Warda, X. 
Vifias, X. Roca-Maza, and M. Centelles, Phys. Rev. C 
80, 024316 (2009). 

C. Xu, B.A. Li, and L.W. Chen, Phys. Rev. C 82, 
054607 (2010). 

L.W. Chen, Phys. Rev. C 83, 044308 (2011). 

Z.Q. Feng and CM. Jin, Phys. Lett. B683, 140 (2010). 

P. Russotto et al, Phys. Lett. B697, 471 (2011). 

E. Chabanat, P. Bonche, P. Haensel, J. Meyer, and R. 
Schaeffer, Nucl. Phys. A627, 710 (1997). 

M. Brack, C. Guet, and H.-B. Hakansson, Phys. Rep. 
123, 275 (1985). 

D. G. RavenhaU, C.J. Pethick, and J.R. Wilson, Phys. 
Rev. Lett. 50, 2066 (1983). 

K. Oyamatsu, Nucl. Phys. A561, 431 (1993). 

J.M. Lattimer and M. Prakash, Science 304, 536 (2004). 

A.W. Steiner, Phys. Rev. C 77, 035805 (2008). 

F. Douchin and P. Haensel, Phys. Lett. B485, 107 
(2000). 

F. Douchin and P. Haensel, A&A 380, 151 (2001). 
J. Carriere, C.J. Horowitz, and J. Piekarewicz, Astro- 
phys. J. 593, 463 (2003). 

S. Kubis, Phys. Rev. C 76, 035801 (2007); Phys. Rev. 
C 70, 065804 (2004). 

A. Worley, P.C Krastev, and B.A. Li, Astrophys. J. 
685, 390 (2008). 

K. Oyamatsu and K. lida, Phys. Rev. C 75, 015801 
(2007). 

C. Ducoin, Ph. Chomaz, and F. Gulminelli, Nucl. Phys. 
A789, 403 (2007). 

P. Chomaz, M. Colonna, and J. Randrup, Phys. Rep. 
389, 263 (2004). 

C. Providencia, L. Brito, S.S. Avancini, D.P. Menezes, 
and P. Chomaz, Phys. Rev. C 73, 025805 (2006). 
C. Ducoin, J. Margueron, and P. Chomaz, Nucl. Phys. 
A809, 30 (2008). 

C. Ducoin, C. Providencia, A.M. Santos, L. Brito, and 
P. Chomaz, Phys. Rev. C 78, 055801 (2008). 
H. Pais, A. Santos, L. Brito, and C. Providencia, Phys. 
Rev. C 82, 025801 (2010). 

C.J. Horowitz and J. Piekarewicz, Phys. Rev. Lett. 86, 
5647 (2001); Phys. Rev. C 64, 062802(R) (2001); Phys. 
Rev. C 66, 055803 (2002). 

H. B. Callen, Thermodynamics, Wiley, New York, 1985. 
J. Xu, L. W. Chen, B.A. Li, and H.R. Ma, Phys. Rev. 
C 79, 035802 (2009); Astrophys. J. 697, 1549 (2009). 
L.W. Chen, CM. Ko, B.A. Li, and J. Xu, Phys. Rev. C 
82, 024321 (2010). 

C. Ducoin, J. Margueron, and C. Providencia, EPL 91, 
32001 (2010). 

C. Ducoin, J. Margueron, C. Providencia, and I. 

Vidafia, Phys. Rev. C 83, 045810 (2011). 

L.W. Chen, Sci. China: Phys. Mech. Astron. 54, suppl. 

I, sl24 (2011) [arXiv:1101.2384] . 

C. Jacoby and S. Nussinov, J. High Energy Phys. 5, 017 



16 



(2007). Astronomical Union. 

[102] Y. Fujii, in Large Scale Structures of the Universe, page [103] K. lida and K. Sato, Astrophys. J. 477, 294 (1997). 
471-477, Eds. J. Audouze et al. (1988), International 



